Measurement of specific contact resistivity using scanning voltage probes 
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Specific contact resistivity measurements have conventionally been heavy in both fabrication and 
simulation/calculation in order to account for complicated geometries and other effects such as 
parasitic resistance. We propose a simpler geometry to deliver current, and the use of a scanning 
voltage probe to sense the potential variation along the sample surface, from which the specific 
contact resistivity can be straightforwardly deduced. We demonstrate an analytical example in 
the case where both materials are thin films. Experimental data with a scanning Kelvin probe 
measurement on graphene from the literature corroborates our model calculation. 



The contact between any two conductors exhibits an 
interface resistance. The intensive quantity that char- 
acterizes this contact resistance is the specific contact 
resistivity, denned as the ratio of the voltage across the 
interface and the current density through the interface. 
The accurate measurement of this quantity is of great 
technical and scientific importance [TH3]. In this Letter 
we propose a method of measuring specific contact resis- 
tivity based on the use of scanning voltage probes. This 
approach is simpler than traditional methods with the 
added benefit that inhomgeneities of the specific contact 
resistivity can be detected. 

The methods that have hitherto been developed in- 
clude, e.g., the transfer length measurement (TLM) and 
the cross bridge Kelvin resistor (CBKR) method. Both 
methods require accurate fabrications of test structures. 
The TLM method 0] extrapolates from two point resis- 
tance measurements of devices with different lengths to 
obtain the contribution of contact resistance in the to- 
tal resistance. Accuracy in fabrication and homogeneity 
are required for this method. In the CBKR method[3HS]j 
a three dimensional test structure is fabricated and the 
contact resistance measured in a four point configuration. 
The measurement result then undergoes a self-consistent 
algorithm to obtain the specific contact resistivity. A ma- 
jor issue when applying these methods is the interplay of 
the length scales of the test structure geometry and the 
materials dependent current transfer length g]. 

In our approach we utilize information from the spa- 
tial variation of the local voltage in the simple structure 
depicted in Fig. [T] to determine the specific contact re- 
sistivity. In this structure there is essentially no geomet- 
ric length scale, and the current transfer length is deter- 
mined from the scanned voltage profile. Both material 
1 and material 2 can be either thin or thick, although a 
thin film might be appropriate for fabrication purposes 
and easier to model. The spatial variation of the voltage 
across the boundary of the two materials can be mea- 
sured, and the data fitted to obtain the value of specific 
contact resistivity. Also, the voltage variation along both 
materials can be used to determine the (sheet) resistance 
of both materials. 



Various instruments can be used to measure the lo- 
cal voltage variation along the sample surface. These 
include scanning potentiometry and scanning electro- 
static voltage probes. For scanning potentiometry, the 
natural choices are scanning tunneling potentiometry 
(STP)[71 [SJ, which has the highest spatial resolution 
(<10nm), and scanning conducting AFM, which is per- 
haps simpler but with courser spatial resolution^ HO] . 
For electrostatic probes, scanning Kelvin probes [TT], and 
scanning single electron transistors [12] (which is a form of 
Kelvin probe with high sensitivity) are suitable choices. 
While both types of probes will work, scanning poten- 
tiometry has the attractive feature that it measures the 
transport potential directly whereas electrostatic probes 
have the complication of dealing with spatial variations 
of the electrostatic potential associated with the patch 
effect. 

The calculational description of the test structure is 
well developed |H [T3]. In our case it is easier due to the 
simpler structure. In the case where both films are thin 
(the condition for a film to be thin will be specified) , and 
sheet resistance and sheet current density can be used 
for both materials, we write down the analytical solution 
below as a demonstration on how this method works. 

The following calculation assumes that the films are 
thin for both material 1 and material 2, described by a 
sheet resistance Rat,2 = Pi, 2/di. 2, where dtp are the in- 
dividual thicknesses of the two films. The specific contact 
resistivity between the two materials is 1Z C . The sheet 
current densities in the two materials are Ji,2, where J2 
is only valid to the right of the boundary [x > 0). As- 
sume that the electrode is macroscopic in size, i.e., the 
other boundary is infinitely far from the boundary near 
which the scan is performed. To the left of the boundary 
(x < 0), due to current conservation the sheet current 
density is a constant: 

Mx) = J x<0, (1) 

and equals the sheet current density applied to the de- 
vice. Thus, to the left of the boundary the local voltage 
is: 

Vi{x) = -RaiJox x < . (2) 
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To the right of the boundary (x > 0), the differential 
equations are: 



dVi{x) 
dx 

dV 2 (x) 
dx 



J\{x)Rui (material 1), 



J-z{x)R\j2 (material 2), 



FIG. 1. Schematic of proposed measurement of specific con- 
tact resistivity between materials 1 and 2. A planar structure 
should be fabricated to present a well defined boundary be- 
tween the two materials. Local voltage variation due to the 
existence of a sheet current density Jo across the boundary 
can be measured and from the measured profile the specific 
contact resistivity can be calculated. A jump in the voltage 
when scanned from above is expected due to the specific con- 
tact resistivity. 
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FIG. 2. Numerical examples of solution j5j. (a) Sheet current 
density redistribution near the boundary between two materi- 
als as calculated in equation |5J. The sheet current density is 
normalized to Jo, the total sheet current density applied; and 
x-direction displacement is normalized to l\ = \/TZ c / Rm- 
(b) Local voltage variation near the boundary between two 
materials as calculated in equation jpj. The voltage V^a;) is 
normalized to JolZ c /h\ and ^-direction displacement is nor- 
malized to h = \/lZ c / Rm- Note that if scanned from the top 
in Fig. [l] the scanning probe picks up Vi at x < but V2 at 
x > 0. Hence a jump in the measured local voltage will be 
observed. 



3\(x) + Ji(x) = Jo (current conservation), 
V 1 (x)-V 2 (x) 



(material interface), 



Jc 



dJi{x) 
dx 



(current conservation at interface), 



with the boundary conditions: 

Vi(0) = 0, ^(0) 



Jo 



(3) 



(4) 



The above equations can be analytically solved. For 
convenience the following two parameters are used: h_2 = 
\JR C / Ru\,2, and we denote I = hh/s/lf + /|- These 
are the transfer lengths in respective materials [4!. The 
solution for x > then reads: 



Ji(x)=J 



J 2 (x) = Jo 
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v 1 (x) = -n c j 



v 2 (x) = -n c j 



72 72 

L 1 L 2 



72/2 

n l 2 



,-x/l 



72/2 

n l 2 



(5) 



Some numerical examples of the above solution, also 
including J\ and V% at x < 0, are presented in Fig. [2] 
The curves show how part of the sheet current density is 
transferred from material 1 to material 2 in the existence 
of contact resistance. The current transfer is the largest 
at the boundary, and exponentially decays with a length 
scale of I. As a result there is a voltage difference between 
the two materials at the boundary, which gradually be- 
comes smaller over a length scale of I as sheet current is 
transferred. 

The solution self-scales, as is evident from both the 
form of the formula and the numerical examples. The 
major features of the solution are: a slope at x < of 
value: 



Pi 



-HcJo/h — —RniJo 



(6) 
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an asymptotic slope at x 3> of value: 

& = -n c j e/{i\il) = -R ai j Q i 2 /il , (7) 

and a jump in the local voltage variation, if measured 
from the top of the device as depicted in Fig. [T] at the 
boundary x = 0, of value: 

Ay = -aKcJo/h = -olBuiJq • h ■ (8) 

Here, a — is a numerical factor that depends on 
the ratio of l\ and I2, as can be seen in the numerical 
examples in Fig. [5Jb). 

The above scaling provides a very convenient pro- 
cedure to determine experimentally all three material 
dependent parameters: Rai, Rn2, an d 1Z C , from one 
scanned local voltage curve across the boundary. 1) From 
the known applied current to the device, Jq can be calcu- 
lated, and Rui = \@i\/Jo can be obtained. 2) From the 

ratio of f3\/p\, R\32 = Rui a r~n a can be obtained. 

\Pl\ ~ \P2\ 

3) From the size of the jump, l\ — |AV|/(a|/9i|) can be 
obtained, from which 1Z C = Railf can be calculated. If 
necessary, curve fitting can also be employed to improve 
quality of calculated results. 

We note the experimental result shown in Fig. 4(d) in 
reference |14j , in which a scanning Kelvin-probe mea- 
surement was performed on a device made of graphene 
(material 1) and Cr/Au (material 2). The authors of 
that paper correctly extracted the contribution of con- 
tact resistance from their data. We would point out, 
however, that this set of data can also be used to estimate 
the specific contact resistivity value. From an estimated 
Ra = 200^/D for graphene and AV/\Pi\ = 2^m from 
their data, we estimate the specific contact resistivity to 
be 1Z C w 1 x 10 — 5 S7 ■ cm 2 , which is very close to values 
obtained in the literature [TS]. 

It can also be seen that I is a characteristic scale in the 
local voltage variation and sheet current density distri- 
bution. Thus it is clear that the film thicknesses should 
be much smaller than I in order for the set of equations 
([3]) and Q to apply. When film thicknesses are larger 
than or comparable to I, a two-dimensional calculation 
is needed to account for current crowding at the edges 
of the electrodes. Nevertheless, the main features - the 
two (asymptotic) slopes and one jump at the boundary 
should be the same. The difference would be a different 
numerical value for a than One exception would 

be the case h ^> h, i-e. a much more conductive upper 
material than lower material. In this long as the 

lower material is thin with d\ <C l\, the upper material is 
an equipotcntial, and the thickness d,2 does not matter. 

To help assess the utility of our approach to the deter- 
mination of the specific contact resistivity, we plot rep- 
resentative values of Rq and TZ C in Fig. [3] The diagonal 
lines in the figure indicate contours of constant current 
transfer length. As seen in the figure, the values of trans- 
fer lengths are between lOOnm and 10/im. This range is 




FIG. 3. Numerical values of 1Z C and /?□ for different materials. 
The lines indicate constant transfer lengths. The markers are 
respective values for: □ NiSi or PtSi to p or n doped Si|16j ; 
■ Si/Ti/TiN to n doped Ge[IZ]; A Ti/Al/Mo/Au to lOOnm 
n doped GaN[TS]; A Ni/Au to GaN nanowireflt]: V Ti/Au to 
ljLtm Al doped ZnO[20]; T Cr/Au to graphene [T4]: Ni to n 
doped 3C-SiCpi): ♦ lOnm Al to lOnm W[22]; and > lOnm 
Cu to lOnm Cu[23]. 



comfortably above the spatial resolution of instruments 
mentioned above. Note also that as the specific contact 
resistivity decreases the high spatial resolution of modern 
scanning probes become more and more essential. 

Practical devices have finite roughnesses ranging from 
nanometers to tens of nanometers [TH [20l [24] . Even 
though the roughnesses quoted here are small compared 
to transfer length values presented in Fig. [3j when the 
roughness is high, there can be a systematic error due 
to the larger contact area than the nominal value, i.e., 
the measured specific contact resistivity is proportionally 
smaller than the true value if the contact region is rough. 
This can be true for both the conventional methods and 
our proposed method. The measured specific contact re- 
sistivity is an effective value not accounting for the rough- 
ness of the interface in both cases. We note that some of 
the proposed instruments in our method, e.g., scanning 
conducting AFM or STP, have the capability of mea- 
suring topographical roughnesses on both materials and 
can thus provide an estimate of the roughnesses, though 
this is not the true roughness of the interface. Systems 
with smooth surfaces and interfaces such as [25l [26] are 
desirable for both contact resistivity measurement and 
fabrication scale down. 

Another finite contribution to the specific contact re- 
sistivity due to sample roughness is the effect of non- 
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uniform current distributions across the interface be- 
tween two materials, across the boundaries between two 
materials, and within the two materials. If the specific 
contact resistivity is small and the transfer length scales 
are comparable to the surface roughnesses, the situation 
is more complicated. These effects will clearly affect 
the accurate determination of specific contact resistiv- 
ity in both the conventional methods and in our method, 
and a simulation is required to address the issues. We 
stress that a scanning voltage probe measurement in our 
method can be used to investigate the deviation from 
ideal curves in the one-dimensional theory, and the two- 
dimensional potential maps on the same sample can pro- 
vide an empirical estimate to the roughnesses. These 
estimates may also help addressing the systematic error 
discussed in the previous paragraph. Such simulations 
required to quantify the effects of roughness will be pur- 
sued in a separate paper. 

Another interesting limiting case is when the sheet 
resistances of both materials are zero R\ji,2 = 0, i.e., 
superconducting. If the sheet current density Jo flows 
through only one material, the case is trivial, namely no 
voltage drop across this material. However, if the sheet 
current is forced into one material and out of the other, a 
voltage drop is expected: superconductivity dictates that 
both materials are equipotentials; because of the current 
transfer, a voltage drop exists at the interface of the two 
materials. If the specific contact resistivity is homoge- 
neous, the current transfer has to be homogeneous over 
the interface, consistent with I — 00. Hence in this case a 
four point measurement across the interface and a good 
measurement of the interface surface is sufficient in de- 
termining the specific contact resistivity 1Z C = AV ■ A/I, 
where A is the area of the interface and I is the current 
through the interface. 

We also note that the problem of thermal boundary 
resistance in the same geometry as Fig. [TJhas a solution 
of the same form with our solution Hence a scanning 
thermometer across the boundary of a thermal interface 
might provide a similarly convenient measurement of spe- 
cific thermal boundary resistivity. 

In summary, the method we propose requires only a 
simple electrode-sample geometry and at the most a light 
numerical simulation. The benefit of self-calibration pro- 
vides high reliability. Multiple instruments are available 
to perform this measurement. 
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